o ; Local Momenta and a Three- Body Gauge 

, Michael Jay Schillaci, Francis Marion University, Florence, SC 29501 

D : 

00 . February 2, 2008 

1 __, ' ' Abstract 

O ,' In recent years researchers have attempted to improve the continuum state three-body wavefunc- 

tion for three, mutually interacting Coulomb particles by including, so called, local momentum effects, 
which depend upon the logarithmic gradient of the continuum, two-body Coulomb waves. Using the 
exact three-body wavefunction in the region where two of the three particles remain close, a revised 
description of these local momenta, is attained and predicts that a quantum-mechanical impulse may 
develop in the reaction zone, causing like-sign-charged particles to decrease their radial separation and 
opposite-sign-charged particles to increase their radial separation. The consequences of these predictions 

• i-h . are investigated through both quantum and semi-classical techniques where the total energy of a two- 

body continuum Coulomb system in the presence of a third, mutually interacting body are analyzed. 
Lh \ Numerical calculations confirm that while ignoring these local effects for light-ion-atom processes, may 

be appropriate, three-body effects may dominate in the reaction zone for heavy-ion-atom processes. The 

i i ' techniques developed here are then applied to explain the observed asymmetry in the data collected by 

Wiese, et. a!.,[PRL 79, 4982] on the correlated breakup of three massive, Coulomb interacting particles. 
The results are attributed to a genuine three-body effect which rearranges particles in the reaction zone, 
while retaining the appropriate asymptotic behavior. Preliminary calculations show a deviation of less 
than 10% between the predicted and observed asymmetries. The strength of these results is then used 
to argue that the local momenta, herein developed, be treated as a formal gauge constraint for three- 
body interactions. This hypothesis is investigated and it is shown that a real- valued, position-dependent 

(*— ■; . phase is added to the wavefunction. A semi-classical analysis of the proposed three-body gauge, reveals 

that while genuine three-body effects may arise in the reaction zone, the asymptotic form of the rele- 
vant two-body Hamiltonian remains unchanged for relative energies greater than ~ leV for all atomic 

^yT ■ species. Further analysis shows however that one may detect asymptotic variations in the scattering 

£_) ' amplitudes for massive systems at energies ~ /ieV. These results provide convincing theoretical and 

" Tji , physical evidence for the success of many current experiments and indicate that more experimentation 

^_ ■ with near-threshold, massive three-body systems is needed. 

>'. 

X: 

Principia) that it was possible to infer the orbits of two, mutually interacting bodies using only the laws 
of mathematics, mankind has endeavored to derive an analytic description of the motion of three, mutually 
interacting particles. 

Three-body interactions abound in natural processes as diverse as stellar evolution and thin-film growth 
and occur over a very wide range of energies. For the present discussion, continuum atomic scattering 
will be considered for three mutually interacting charged particles interacting via the infinite range Coulomb 
potential. Calculation of local distortion effects herein derived will be carried out over representative energies 
of between 3 - 54.4 cV.Q 

1 These energies were chosen because of the particular relevance to the results offered in [pj on electron or positron scattering 
from hydrogen. 
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1 Introduction 

The "three-body" problem is as old as the study of Physics itself. After Sir Isaac Newton showed (in his 



In Section 0, the traditional Jacobi coordinates are introduced and key historical results are given. 
While the notation used here is not substantially different, minor errors have been corrected and additional 
properties and unique results have been added. 

A revised description of the so called local momenta, first presented in [pi is offered in Section 0. While 
the local momenta derived here also depend upon the logarithmic gradient of a continuum state, Coulomb 
wave - here referred to as the local distortion - the coordinate dependence is such that variations of the 
momenta can not be ignored in an a priori manner. The development continues in Section H, where a 
detailed discussion of the mathematical and physical properties of the local momenta is offered. Appendix 
|A| shows that the local distortion can be expressed as a damped oscillator function that is analytic over all 
regions, thus improving the ability to assess the possible contributions of these local effects in both the inner 
(reaction-zone), and the outer (asymptotic) regions of the three-body scattering event, [j This analysis reveals 
that while local distortion effects generally "fall-off" in the asymptotic regions, they may alter the outcome 
of a scattering event by rearranging particles in the reaction-zone. Appendix || demonstrates however that 
one may incorporate local effects while retaining the appropriate asymptotic form for the continuum state, 
three-body wavefunction. From this a new interpretation of the three-body scattering event emerges, wherein 
a two-body continuum pair acquires a local momentum by scattering off of an exact interaction potential. 

Section g illustrates that generalized local momentum effects may be used to explain the observed asym- 
metry in the data collected in the triple-coincidence detection experiment of Wiese et. al., [|2|. The experiment 
considered three large, very nearly equal mass, mutually interacting, charged particles. With the appropriate 
total center-of-mass energy and reduced mass parameters chosen to reflect those used in the actual experi- 
ment, the predicted local effects are shown to be large enough to reproduce the observed asymmetry. This 
result is taken as motivation to hypothesize that the generalized local momentum be treated as a possible 
gauge transformation for three-body interactions. In Sec. o this hypothesis is formalized and a semi-classical 
analysis reveals that while a real- valued, position dependent phase is added to the two-body continuum state 
wavefunction, the asymptotic form of the relevant two-body Hamiltonian remains unchanged. 

2 Notation 

The traditional Jacobi coordinates, here denoted by (r„, p v ) along with their respective, congugate momenta, 
(k v ,q v ), where v = a,/3, 7 will be used to indicate a particular channel representation. These coordinates 
are particularly well suited for the study of the motions of three mutually interacting particles because they 
locate the conventional] reduced masses of the system. These coordinates are shown in Figure [j] and the 
"alpha-channel" representation will be used unless otherwise specified. 
The reduced masses located by the (f a , p a ) coordinates are given by 

M« = ^- J - (la) 

77i/3 + rn 7 

M a = "«("/» + "7) (lb) 

m a + nip + m~/ 

respectively and the other masses are defined cyclically. 

In addition, there are two relationships that the Jacobi coordinates and their conjugate momenta obey 
that will be of use in the current development. These are, 

r a + f/3 + f 7 = (2a) 

ho.. + hit + h. = 0. (2b) 

2 While there is no rigorous convention for the quantification of these regions, the "reaction zone" is here taken to be the 
region wherein r < 10ao. 

3 While there is some debate in the literature, as noted in |H|, over the appropriate reduced masses to use in the Classical 
treatment of the three-body problem, the possible renormalization of mass in the Quantum treatment makes these arguments 
irrelevant here. 



While (2a) can be "seen" in Figure (jl, (Eq) is more subtle and is a statement of the "relative- velocity 
conservation" for the three-body system. These relationships follow from the orthogonality of the Jacobi 
coordinates and are easily verified using the transformation matrices for the coordinates, 
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and for the momenta, 
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A further consequence of these transformations is that the three-body, Coulomb potential may be ex- 
pressed as follows^: 
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where the Z v [y = a,/3,j) are the appropriate charge signs. Moreover, the orthogonality of the Jacobi 
coordinates can be used to show that, the three-body Schrodinger equation may be written in a channel- 
independent way. i.e., 



(HZ + VS) %^ (?„,&) = Et*^ (fv, P.), 
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and Vq(t v , p v ) has been used, with supressed coordinate dependence, to indicate the full three-body Coulomb 
potential of the z/-channcl. Here £ t is the total center of mass energy of three-body system. 

To date, the most successful and widely used approximate solution to the three-body Schrodinger equation 
is the paradigm "3C" wavefunction proposed by Redmond @[,[g| and rigorously derived and tested by Brauner 
et. a/.H. The solution is valid in the asymptotic region where all particles are far apart. Traditionally one 
denotes this region with, Qq, and the solution is given by, 
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where for instance, 
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4 Whilo it is conventional to use "atomic units" wherein, ft = e = 1, all physical constants are retained so that the interested 
reader may verify explicit numerical results cited later in the text without normalization. 



is the atomic Sommcrfcld parameter. The hyper-parabolic coordinate, 

Cfc„ = k v r v + k v -rv, (10) 

has been introduced for notational simplicity only and with the above definition, the wavefunction (H) 
satisfies all incoming boundary conditions in the region Qq. The overwhelming opinion in the literature is 
that all valid three-body wavefunctions must match smoothly with this solution in the region £~Iq. 

The logarithmic phase factors in the Redmond solution are present physically because of the infinite 
nature of the Coulomb potential, and arise mathematically as the leading term in the asymptotic expansion 
of the confluent hypergeometric function, i.e., [|7| 

<V(CO = W„^i*i(«7fc„l,-*CkJ (Ha) 

C%„(CO = r[l-i7y fc Je 2 2^— : (lib) 

n=0 ( L )n n. 

C m (CO = I™ e"""- 1 "^. (lie) 

r u — >oo 

These "C-functions" provide one of the representations for the exact solution to the two-body scattering 
problem. That is, if one writes the two-body wavefunction as 

^°- Bodv {r l/ ) = e^CV„ (CO, (12) 

then substitution of this form into the two-body Schrodingcr equation shows that 
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While the very intuitive solution m) has been used with great success by researchers to model both atom- 
ion m,M and photo- ionization processes |fi0[ , a more robust solution has been sought in recent years M,[E|. 
Particularly, a solution is sought that may be extended into the so called "interior regions," where at least 
two of the three particles remain close to one another. It is significant to note that all of the successes of 
the 3C wavefunction have been achieved for light-ion - atom systems, where the asymptotic form in Qq has 
been shown to be generally adequate. What is sought however is a more precise accounting of the intimacies 
of the three-body interaction for arbitrary masses and in all regions of the scattering space. 

3 Origin of the Local Momenta 

Because an exact solution to the atomic three-body problem does not exist, the best hope in achieving an 
improved wavefunction has been to improve the approximation schemes used. Generally, these approximation 
schemes fall into two categories: 

• Approximating the Kinetic Terms of the Hamiltonian (i.e., the Eikinol approximation) 

• Approximating the Potential Terms of the Hamiltonian 

While these techniques would seem to be mutually exclusive, the problem is that the inseparable nature of 
the three-body Coulomb potential makes the choice of the kinematic description impossible. Research has 
thus continued along a fragile path and to account for approximation and/or distortion effects, two well 
established interpretations have emerged: 

• Introduce a Local Momentum, which depends upon the radial separation of two of the three particles 
through the logarithmic gradient of the continuum two-body wavefunction, and attribute the distortion 
to a velocity-dependent, auxilary potential, cf., O 



• Introduce an Effective Charge, which depends upon the radial separation of two of the three particles 
through the logarithmic gradient of the continuum two-body wavefunction, and attribute the distortion 
to dynamical screening of charges, cf., pi 

The common dependence upon the two-body solution in these two interpretations is clear and in both 
the dependence is derived in a rigorous, but a posteriori way, to satisfy the relevant boundary conditions 
in the asymptotic regions. The natural question is if this dependence can be achieved in an a priori way, 
utilizing only physical and mathematical intuition. 

To approach an answer to this question, it is important to know that there does exist an exact solution 
to the continuum three-body Schrbdinger equation for the case when two of the particles remain close 
together(or equivalently if one of the particles is infinitely massive). The solution depends critically upon 
the form of the three-body Coulomb potential in the region fi Q , where two of the particles, (3 and 7 are close 
together, while particle a is infinitely far away from both of them. That is, one may write a description of 
the region f2 Q succinctly as follows: 
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Clearly the region f2 Q matches smoothly with the region fi , and to investigate the behavior of the three-body 
Coulomb potential in this region, one may write 
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where 
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and the Ph'^Ja • pa) are the Legendre polynomials of the first kind. 

Now in the region O a , the only surviving term in the expansion is for V = 0. Hence the potential has 
the following asymptotic form: 
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The second term in this expansion is often referred to as the "reduced charge potential," and the resulting 
form of the three-body wavefunction is separable, e.g., the three-body Schrodingcr equation takes the 
asymptotic form, 
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separates and if one assumes that the three-body wavefunction has the standard plane wave form 
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where the total energy of the three-body system may be written as 
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for continuum scattering, then substitution into the asymptotic three-body Schrodingcr equation yields the 

Z^Zp + Z^e 2 1 



following equation for $£ - (f a ,p a ) in the region f2 c 
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Therefore the exact solution in the region Cl a is given by, 

*£1(W«) = e^^+^C^ (Ck a )C Vqa (C S J. 
Here the function C,, (£ 9q ) is a reduced charge continuum state Coulomb wave, which satisfies 
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While equation ( J20| ) constitutes a rigorous solution it is not a valid solution because it does not match 
smoothly with the result asserted by Redmond, (cf., equation (H) above.) To sec this note that the asymptotic 
form of (pOl) in O would be given by 
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Clearly, the two logarithmic phases present in this solution can not match smoothly with the three present 
in the Redmond solution. Hence ( |20|) does not constitute a valid solution in the asymptotic region f2o- 
In addition, note that as a general wavefunction in the region fig. (|20| ) may not be a good choice simply 
because when the two charges Zp and Z 7 are of equal and opposite charge, as in the paradigm case of 
electron scattering from hydrogen, the solution in (T a reduces to that of a free particle, i.e., 
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These reasons make it clear that in order to extend the 3C wavefunction into the the asymptotic region, 
£l a (or into the interior regions for that matter) a different approach is warranted. The current approach, 
refered to here as "kinematic coupling," involves the introduction of an exact interaction potential given by: 
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Hence the three-body Schrodinger equation takes the exact form 
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To remain completely general, one then asserts that the total center of mass energy may be written in 
the form 



£t = £a a + £m + £i , 



(26) 



where £ fla = 2 ° > and £u a = -^f- > are the energies associated with the continuum state, two-body 
clusters, p a and M a respectively, and the term £j accounts for the remaining energy of the three-body 
interaction. Though the definition ( |26| ) is nonstandard, it reflects the fact that because the three-body 
potential must remain inseparable in a completely general solution, so too must the total energy. 

Due to the fact that the inseparable portion of the three-body Coulomb potential is now contained within 
the interaction potential, one can assume that the three-body wavefunction takes the form, 
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(27) 



where Xu - (r a ,Pa) is an unknown function, and C Vk (CO an d C,, (Cg Q ) are the known Coulomb waves 
defined in equations dl3) and (|2l|) respectively. That is, the ansatz (E7J) incorporates the exact solution in 



the region Q, a (cf., equation (2C)) explicitly. 

Substitution of the form (27) into the three-body Schrodinger equation (J25J) yields the following result: 
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where the exponential terms have been cancelled from both sides of the equation and the following vector 
identities have been employed: 

V(AB) = BVA + AVB, 

and 

V 2 (AB) = BV 2 A + AV 2 B + 2VA ■ VB. 

Using equations (|l3| ) and (|2l|) in conjunction with the definition (pq), one finds that the first two lines 
in this unwieldy expression are identically zero. Then after dividing by the product C Vk (Ck a )^vq a (Cg Q )i the 
following equation for the unknown function, x% - (fa. Pa) is derived: 
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is found, where 
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are the proposed position-dependent local momenta for this kinematic coupling model. 

The principle difference in this development is that the local momenta arise purely from the physical 
structure of the inseparable three-body Schrodinger equation, and that their mathematical form implies the 
existence of position dependent momenta. Furthermore, in the current development, the local momenta are 
predicted in a symmetric fashion so that both "legs" of the three-body interaction experience distortions 
which depend upon the congugate coordinates. That is, f a is congugate to k a and p a is congugate to q a . This 
means that the distortion that two of the particles experience does not depend explicitly upon the distance 
to the third particle as in B , but is rather wholly attributable to a genuine three-body effect wherein a very 
intuitive description of the three-body scattering event arises: 

The two reduced mass clusters p a and M a , initially described by the two-body waves, C Vk (Cfc Q ) 
and C' v (C<j„) respectively, scatter off of the interaction potential Vi(f a ,p a ) and acquire a local 
momenta. 

The subsequent motion of these clusters is then of course dictated by the function \% - (?«) Pa) which must 
satisfy (P9|) , 

While one may argue that solving equation (g9|) is a greater task than solving the original Schhrodingcr 
equation, a solution that is consistent with the 3C wavefunction in the region fig may be established. The 
details of this solution are not important to the present discussion and are relegated to the appendices. (See 



appendix fBu What will be of great importance is the nature of the predicted local momenta, (30) 



4 A Generalized Local Momentum 



The local momenta introduced in this kinematically coupled model depend upon the congugate coordinate. 
Because of this variations in the momenta will generally contribute to the solutions. To understand how and 
where these contributions will be important, a generalized "local distortion" term, Dp{f) is introduced, where 
p and r may be either k v or <f„ and r v or p„, respectively. Hence a general position-dependent momentum 
may be defined as follows: 



P(r)=p-Dp(r), 
where the exact form of the local distortion is given byn: 
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and (5,j p (Cp) is a reaZ-valued, position dependent phase, defined by 
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While the general form (32a) has been used by many researchers, the form (32b) is unique and has been 
introduced to help illucidate the physical significance of the local momentum, e.g., one may write: (r ^ 0) 



Dp(f) = Z 



M a n P 



(C P ) 



(p + f), 



l-Kf-ohJ V 2S 

where Z and \x are the relevant product charge and reduced mass of the pair described by C v (C, p ), and 
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The form ( pq ) shows that the local momentum is essentially a - damped-oscillator function, and that 
it is analytic everywhere except possibly at r equal to zero. Interestingly, it can also be seen that the local 
momentum has equal radiation (in the f direction) and induction (in the p direction) components. This is 
indicative of the possibility of a tensor force in a Classical treatment or of a nonconserved current density 
in a Quantum treatment of the three-body interaction. Perhaps more importantly, one sees that the local 
distortion depends upon not only the radial separation of the constituent particles, but also upon both the 
relative energy and reduced mass of the system. 

In Figure the real part of the distortion experienced by a electron-proton continuum pair has been 
plotted with various relative energies and an arbitrary scattering angle of 9 = 0. One sees immediately that 

5 See Appendix IN for details concerning the derivation of this form. 



the distortion "falls off" very quickly with increasing radial separation and even more dramatically with 
increased relative energy. 

To see how the distortion depends upon the reduced mass of the system, observe that in Figure pi the real 
part of the distortion experienced by an electron-electron pair and a electron-proton pair have been plotted. 
Notice that the distortion increases accordingly with increased mass, and that the sign of the distortion is 
different for the cases of attraction and repulsion. 

This very interesting point will be discussed in greater detail below, but for now note that in addition to 
the important physical properties of the local distortion, there arc also important mathematical properties. 
The most important of these is that because the local distortion is essentially an oscillatory function, the 
existence of turning and stationary (i.e., maxima and minima) points may yield much insight. 

In Figure |j one sees that the stationary points^] occur quite often, so that the local distortion will very 
often be identically zero! Moreover, these zeroes will be very dense on a macroscopic scale and so may 
greatly alter the topology of a scattering event. This can be seen even more dramatically in Figure |j[ where 
three-dimensional and contour plots of the local distortion have been shown. 

Though there is no analytic form for determining the zeroes of the confluent hypergeometric function, 
the stationary points of the local distortion may be found by finding the nonzero roots to (|3J). Hence one 
ends up solving the following transcendental equation: 

l = w^m {1 - cos9) > (38) 

where n is an integer. Note that if S p (f) were a constant, (pq) would be the equation of a conic section with 
an eccentricity of 1 - parabolas. Hence one sees that the stationary points of the local distortion lie (very 
nearly) along the classically forbidden "trajectories" for particles with positive energy! While this may seem 
to be no more than coincidence, 5 p (r) is in fact very nearly constant between the stationary points (see Fig. 
0), and one may infer that these "zero-distortion-trajectories" are in fact the dominant contributions in a 
path-integral approach to the quantum, three-body scattering problem. Indeed, from this point of view, one 
may argue that they must be the dominant contributions for light-ion-atom processes in the asymptotic 
regions, where great success has been achieved while ignoring local distortion effects. What is not clear 
however is how these local effects will contribute in the so called "interior regions" or for heavy-ion-atom 
processes. 

5 Kinematic Rearrangement 

Above it was noted that the local distortion had different signs for the attractive and repulsive cases. The 
real importance of this observation is that if one imagines a three-body system composed of a tightly bound 
continuum state pair, with an initial momentum p, and a third, mutually interacting particle, then upon 
break-up the continuum pair will acquire a local momentum P(r). Hence the continuum pair will experience 
a momentum change (or impulse) given by 



1(f) = $tP(f) - p = -$tD f (r). (39) 



Evidently one finds, (for r ^ 0) 
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The equations ( f40| ) imply that the two particles that are initially attracted to one another, within the 
continuum pair described by C Vp (Cp)i wu l experience an impulse that will tend to increase their radial 

6 That is, the points where the zeroes of the real and imaginary parts of the local distortion coincide. 



separation upon breakup and conversely, that two particles that are initially repelled from one another will 
experience an impulse that will tend to decrease their radial separation. In other words, at small values of 
the radial separation, opposite(same) sign-charged particles will tend to repcl(attract) one another due to 
the local distortion effects! 

To illustrate this point, the real part of the distortion for several different continuum pairs have been 
calculated and the results are shown in Figure u\. Observe that even though all of these distortion terms 
decrease with increasing radial separation, the magnitude of the distortion in the interior regions depends 
critically upon the system being studied. As was shown above in Figure 0, the distortion effects increases 
dramatically with an increase in the reduced mass. Specifically, note that for the case of a proton-anti-proton 
((e) if Fig. 0) or proton-proton ((f) in Fig. 0), the range has been shortened accordingly, to show that the 
distortion in the reaction zone is dramatically changed. Indeed, for the case of the proton-anti-proton, 
the magnitude of the local distortion effect may be large enough and in a direction such that a kinematic 
rearrangement of the particles may occur. In addition, because of the dependence upon the inverse of the 
square-root of the relative energy, these effects will be even more pronounced at lower relative energies. 

This kind of phenomena can be observed in the data obtained by Wiese et. al., [|2| on the breakup of the 
excited ion (H^T )*. This highly unstable ion of hydrogen decays in a two step process, as follows: 

(H+y — > H%* + H+ -> Hf + H-+ H+. 

The data (see Figure g) show a slight asymmetry with small increases in energy (from 6.5 eV to 7.5 eV) 
which can be accounted for in a direct manner by considering the local distortion effects herein derived. 
To see this note that during the nearly co-linear breakup of the excited H^f ion the "initial"^ proton Hf 
is emitted in the forward direction and the doubly excited H%* ion is emitted in the backward direction. 
The subsequent breakup of this ion occurs such that a "final" proton Hf and the nearly equally weighted 
H~ ion are formed. M At this point the H~ is preferentially associated with H f . However, after the 

breakup, the continuum pair (H~f,H~) will acquire a local momentum, P(r), which has a component in 
the r direction. Hence this gained momentum will act to increase the radial separation and will push the 
H ~ over the Coulomb Saddle, so that it H~ will be preferentially associated with the "initial" proton, Hf. 
Now, because the two protons are in fact indistinguishable, one will observe an asymmetry between "low" 
and "high" energy scattering events. To help visualize this process, a schematic diagram (sec Fig. pi) has 
been developed. 

Further analysis shows that, as mentioned above, because the local distortion depends upon -4=, as the 
energy is increased, a proportionately smaller number of particles will experience an impulse that is large 
enough to alter their final state distribution. Indeed, during the breakup of the Hf ion, the magnitude of 
the local distortion effects experienced by the (Hf , H~~) continuum pair while in the reaction zone arc many 
orders of magnitude larger than the effects experienced by an electron-proton continuum pair, (see Fig. nfl) 
With this understanding, one may conclude in a straightforward manner, that the observed asymmetry can 
be attributed to local distortion effects of the three-body system. Moreover, the degree of asymmetry may 
be predicted as follows: 



-Ehigh amplitude of rearrangment at low energy 

— ex . (41) 

-Blow amplitude of rearrangement at high energy 

Therefore, with the probabilistic interpretation of the wavefunction, one may write, 



f^high _ / ' low ,.n\ 

-^low y ^nigh 

where 7\igh and Plow are the probabilities of rearrangement for the high and low energy states, respectively. 
That is, using the actual data (See Fig. g.) the ratio of the number of rearranged particles to the total 

7 "Initial" and "final" are used here as in [G to indicate the order in which the particles were detected during triple coincidence 
measurements. 
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number of detected particles may be calculated. Doing this one obtains the following results: 



ghigh = / 7.5eV = 
Slow V 6.5 eV 

and 

/n 



1.2 



V 47 

The absolute error between these two calculations is ss 8.7% and illustrates that the generalized local 
momentum, herein developed may in fact be the actual momentum of the (H~ , Ht) continuum pair during 
breakup. Indeed, one may consider revaluating the importance of the generalized local momenta. Instead of 
treating it as a mere mathematical nicety, one may place it on firm physical ground by viewing it is a formal 
gauge condition for three-body interactions. 

6 Towards a Three-Body Gauge 

If one grants that the mathematical and computational evidence gathered here, is sufficient to guide further 
experimental and theoretical investigations by considering the generalized local momentum as a formal 
three-body gauge condition, then one can construct a three-body gauge transformation, i.e., one may write: 

P^ P(r)=p-Dp(f), (43) 

when working with three-body systems, in much the same way that one would write, 

p — > P{r) = p- ehA(f), (44) 

when working with particles in an external magnetic field, (i.e., the Coulomb Gauge.) 
The subtlety here is that while the transformation, 

A(f) — ► 1(f) + Vfft(r) 

preserves the form of the magnetic field, B(r) = Vp x A(r), and provides for the global gauge symmetry of 
the electromagnetic interaction, the condition (E3I) will contribute a reaZ-valucd, position-dependent phase 
to the relevant wavcfunction.n Thus any proposed three-body interaction mechanism would require a local 
gauge symmetry.fj 

To see how a position-dependent phase arises, one may consider a situation similar to that discussed in 
Sec. |5|. That is, consider three-body system consisting of a continuum pair described by the wavefunction, 
^p(^) = e in ^ r ^vp{^,p) together with a third, mutually interacting particle. Then, upon breakup, the pair 
will acquire a local momentum and according to (|43|), the wavefunction will transform as follows: 

%(0 — *P&<?> = ^ tlP{rl ' rl C Vp (C P ) = e^e+* 3 ^>r CvpiCp y (45) 

Using (ph] ), one then finds thatQ (for r ^ 0) 

Dp{r) ■ f = a T] p [l - cos(C P + 25 Vp )] + ia r) p sin(£ P + 26 Vp ). (46) 



8 It is interesting to note however that, Vp X Dp(r) = «Vp X VplnC^ p (f p ) = 0. 

9 For those readers that are unfamiliar with the terms global and local in this context, note that a global gauge transformation 
is one that preserves the modulus of the appropriate wavefunction. eg., \P(f) — > e~ m \&(r). On the other hand, a local gauge 
transformation generally does not and may take the form, \t(r) — > e^(^\I/(r). 

10 Note that the "extra" ao arises here because of the scale used in the plots, i.e., one writes r — + aorr. 
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Therefore the transformed continuum two-body wavefunction takes the following form: 

\jj ^ ff\ — e ip-r e ia, r] P [l-cos(( p +2S rlp )] e -a a 7 lp sin(( p +2Sr, p )(J f{\ f£j\ 

which shows the explicit form for the position dependent phase, i.e., one may write, 

S n P (C P ) = ia-oVp t 1 - cos(C p + 25 Vp )] - a r] p sin(C p + 25 Vp ), (48) 

and see that the transformation is indeed indicative of a global gauge symmetry, eg., 

^ (f) (f) = e s "^«%(r) (49) 

It is of course the absolute square of the wavefunction that is truly important for predicting whether or not 
this phase will significantly effect the experimental findings. To this end, one may take the absolute square of 
(E7j) and show that the imaginary-part of the local distortion contributes a reaZ-valued, position-dependent 
phase, i.e., 

|^(r)| 2 = e -^on P s in (c P +2S Vp ) \c Vp (( p )\ 2 = e 2 ***^ \C Vp (( p )\ 2 . (50) 

The exceedingly small magnitude of the term 2aoi] p sm (Cp + 2S Vp ) for the relative energy and reduced masses 
of the systems considered in current research findings, makes it clear that one may write 

e -2a J)pSin(Cp+2(5, p ) > ^ /g-j\ 

While this reinforces the fact that the proposed framework will leave the asymptotic description of the three- 
body scattering event unaltered, (as was established in Appendix |BJ) it docs not address the extent to which 
the proposed gauge transformation, (|43|) will alter the description in the reaction zone. 

To begin an investigation of the expected behavior in the reaction zone, one may construct a semi-classical 
expectation value for the total energy of the continuum pair, C rip (C P )- To do this recall that the expectation 
value of the semi-classical Hamiltonian for this system (before breakup) would be, 

^'^-^r+i^r- (52) 

where p = ^ is the momentum of the continuum pair before breakup. During breakup, the pair will 
acquire a local momentum, so that the expectation value of the semi-classical Hamiltonian for the system 
would then become, 

The effects of the proposed gauge transformation may then be observed by plotting the expectation 
values and varying the relative energy, £ (see Fig. |ll|), the scattering angle, 9 (see Fig. [L2|), and the reduced 
mass, [i (see Fig. [F|). All of these show undeniably that while the asymptotic form remains unchanged 
with respect to variations of all kinematic parameters, genuine three-body distortion effects may arise in 
the reaction zone. Note specifically that the variation of the relative energy with the scattering reinforces 
the interpretation offered in Sec. ||. Specifically, the magnitude of the distortion effects for large scattering 
angles (corresponding to the nearly colinear breakup of the continuum pair) in the reaction zone may cause 
two oppositc-sign-charged particles to increase their radial separation, and appear to repel one another. 
Conversely, two like-sign-charged particles would decrease their radial separation, and appear to attract one 
another. 

As a last exercise, one may ask if the position-dependent phase introduced by the three-body gauge 
formalism could be measured. To answer this question one may observe that 



V2z, 



e -lr] p aosm(Q p +16 p ) ^ e -M lp a _ g - ^ lQh y ^g ^N 
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where (p3J) has been used for r] p , and £ is measured in electron-volts. As mentioned above, this contribution 
approaches unity for all systems of physical interest in the realm of current experiments in atomic physics. 
One can however venture outside of this realm and ask what energy and/or reduced mass is needed to obtain 
a measurable deviation from unity. If one uses the heaviest purely atomic species, cither a proton-proton 
ip~ p) or a proton-anti-proton (p — p) continuum pair, and assumes that a deviation of one part in a million 
can be measured, then the energy scale needed is on the order of 1 /icV! The variation of the term e _2aor?p 
at these energies is shown in Figure [14, and illustrates that one may detect a change in the absolute square 
of the wavefucntion, thus altering the relevant scattering amplitude. Moreover, the exceedingly small energy 
scale required to detect these asymptotic distortion effects in either electron(or positron) scattering from 
hydrogen or in electron-electron or electron-positron ionization processes provides a precise understanding 
for the success achieved in these areas while ignoring local momentum effects. 

7 Conclusion 

In the above it has been shown that a generalized, position-dependent local momenta, which depends upon 
the congugate coordinate through the logarithmic gradient of a continuum state Coulomb wave, may provide 
evidence for the manifestation of genuine three-body distortion effects in the reaction zone. The form of this 
local momenta was derived from a consideration of the exact three-body wavefunction, for the case when 
two of the three, mutually interacting, particles are far apart, and it indicates that the effects do not depend 
explicitly upon the location of the third particle. For this reason, the effects may be viewed as a distortion 
of the initial two-body continuum state wavefunction of the two remaining particles. This interpretation was 
adopted, and it was shown that the local distortion effects could be used to provided a rigorous, physical 
description of the observed asymmetry in the data obtained by Wiese et. al., Q on the breakup of three 
massive Coulomb particles. The degree of this asymmetry was then predicted with an error of less than 
10%, and it was also shown that while the distortion effects were large in the reaction zone, the asymptotic 
form of the relevant two-body interaction may be retained by treating the local momentum acquisition as 
a three-body gauge constraint. Furthermore, the evidence for detecting asymptotic variations in the /xeV 
range, as presented in Fig. n4 suggest that more experimentation be focused on these low energy, heavy- 
ion-atom processes. Indeed, these experiments may yield new insight into a mechanism by which electrical 
forces may contribute to the fusion process! That is, the quantum-mechanical-impulsc interpretation offered 
here shows that two like-charged particles may in fact attract one another due to local distortion effects in 
the reaction zone. 

In addition to these findings and predictions, one may learn much by noting that by adopting the proposed 
three-body gauge transformation, one finds that an electron-proton continuum pair exhibited a very small 
amount of distortion in the reaction zone. This result provides a rigorous explanation of why the paradigm 
3C wavefunction works so well for light-ion-atom scattering ra,B. Moreover, the general framework shows 
that an electron-electron continuum pair would experience distortion effects of lesser magnitude, due to its 
greatly decreased reduced mass. These results can again be used to explain why Qiu et. al.,Q achieved 
amazing success in modeling electron-electron photo-ionization processes, while ignoring local momentum 
effects. Indeed, if one recalls the path-integral interpretation suggested in Sec. |5| then the success of these 
findings for light-ion-atom processes may be attributed to the fact that the leading contribution to the 
relevant cross sections are the "paths" along which the local distortion effects are identically zero. While no 
analytic form exists for calculation of these roots, one can compile a table for use in numerical calculations 
and construct a solution that better reflects the physical nature of the three-body interaction. 
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A The Logarithmic Derivative of the iFi(if] p , 1; — l£ p ) Function 

The logarithmic derivative of the confluent hypcrgcomctric function is defined here as follows: 

D p {f) = i$ r ln 1 F 1 (ir) p ,l;-i(p)=i rl ]. — P ' ' f 

(55) 

iFi{ir] p ,l;-i( p ) 
Then using the well known recursion relationship, (See reference M for instance.) 



fi.Fi [{a + 1), (b + 1), z] = iFj [(a + 1), b, z] - X F X (a, 6, z) , (56) 

o 



one may write] 1 j 



D P <r) = iW [lFl ^ + f % ~ iQ 7 ^ 1; - %Cp)] (P + f ), C P >0. (57) 

Further simplification of the logarithmic derivative is obtained by use of the Kummer relationH for the 
confluent hypergeometric equation, e.g., 

!F! (a, b, z) = e\F 1 [{b - a), b, -z] . (58) 

With this relation one finds that the distortion may be rewritten as follows: 

f> ,~ ■ [e- lQp iFi{-ir 1 p,l;i^ p )- 1 F 1 {iT lp ,V,~iC P )\ ^ 

D P {r) = ir) p p — — ■ (p + r) (59a) 

Cp iFi(ir] p ,l;-iC P ) 



r p ■ (p + r) 



(P + r), r > 0. (59b) 



\ _ c - t c„ i F i(- ii7 ?P' 1 ^Cp) " 
1 F 1 (ir] p ,l;-i( p )_ 

11 Because of the division by <^ p in equation (fe7J), the region of validity for thclogarithmic derivative is limited to £ p > only 
as indicated. One can however find the actual value at ( p = 0. c.f. equation (B3). 
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At this point note that the ratio of confluent hypcrgcomctric equations in equation (59b) is in fact a very 
special case because the function on top is the complex conjugate of the function on the bottom! i.e., 



i-Fi(-iT?p,l;«'Cp) = iF^(in P ,l;-iC P ), 
so that one may define a reaZ-valued, position-dependent phase, 

tan \S (C )] - Sti^ifofr. 1 ;-^)] 
The logarithmic derivative may then be written as follows: 



Dp{r) 



iVpPiP + r), r = Q 



>ip 



I _ e -*(C P +2<5„ p (Cp)) 



(P+f) 
p-{p+f) ' 



r >0 



(60) 



(61) 



(62) 



and one sees that in this form, the local distortion is analytic everywhere, except possibly at r = 0. This 
result is very significant, because it allows one to investigate the nature of the local momenta well inside 
of the interior regions of a scattering event in a rigorous manner. Only in this way can one determine the 
relative importance of these local effects. 

B An Alternative Three-Body Wavefunction 

To solve equation ( |28|) f or the unknown function \t - {r a ,Pa) one normally assumes that the interaction 
energy (cf, equation (Pq)) satisfies 

£/ = 0. 

This technique was developed by Popalilios as referenced in 0, and asserts that the total center of mass 
energy of the three-body system is partitioned among the two reduced mass clusters, fi a and M a . Hence 
one may write 



■JL-V 2 Pa -i— K(r a )-V Pa 



ft h 



XL„-(' ; 'a,P«)=0. (63) 



and see that Xj: - (fa. Pa) must be built to incorporate each term in the interaction potential. 
To accomplish this, assume that Xk - ('qi Pa) is given by 



Xfc g {r a ,p a ) = C^ lqa {C qa )fp{r cn p a )f 1 (r a ,p a ), 



where C- Vqa (( qa ) satisfies 



" V p-o " *TT& ' V p» 



Z a (Zj3 + Z 7 )e 2 1 



c Vq (C q J = o. 



(64) 



(65) 



At) 



Because the function C-^ (C, qa ) is a known function, which exactly incorporates the interaction term 

A-KtoPa 



- — , this substitution results in the following (exact) coupled equations for the functions f v {r v , p v ):(y 



<L V l ~ Yh~y% ~ ijtK2(?a;?y,P-y) • V r -„ - i£Q2(#„ *,,/*,) ■ V p -„ + V™ a (fp,p a )\ W ,pp) = 0, 



h' V72 _ h' V72 
2(i„ v r„ 2A/„ 



I'.vn 

:*}« ~ ^K^(f a -f ,p fi ) ■ Vr, - ifcQgfa,?/,,?!,) ■ V p - a + V^%(f 7 ,p a )] / 7 (r 7! p 7 ) = 

where (for e^(/, with e = j3, 7) 

K^(r a ;r e , p e ) = k a - iV? haC %Q (OtJ - -Vf a ln/ £ (r e , p e ), 



(66a) 
(66b) 

(67a) 



Q%(p a ; r £ , p £ ) = q a - iXlp a [lnC, qa (Cg a )C_^ a (C 9 J] - ^/? a ln /e(r e , p e ), (67b) 
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idQ 



2 



,&(,„„_ «£± X { 1 + J!!3!S 



47reo r u Z a Z e e 2 M a \p a 



r„ 



1 _ e -KC qa +2<5,„ jQ (C,„ )) i_ e -i(C qa +2-5-„„ a (C,« )) 



(68) 



While these equations may seem even more complicated than the original equation, there are many 
important aspects to note: 

• The effective potential reduces to the Coulomb potential for the ^-channel, in the asymptotic regions 
fio and Q a - i.e., 

v (e)(n ,n a ) ,-. -, s __ Z a Z e e 2 1 

• The second term in Q £ a {p a ,r ei p e ) vanishes in the asymptotic regions because, 

[C Vqa (C 9 JC-„ ta (C,»)] !!»,„_ — e-*»«- 1 "<-e+"'- ln ^ = 1. (70) 

(Indeed, this cancellation of the logarithmic phases was "built in" to the solution!) 

• The first two terms in K^(f a ,r £ , p £ ) may be identified as a position-dependent, local momentum, i.e., 

K a (r a ) = k a - iVp a h\C Vka «*„,)• 

• The solutions are coupled in a completely symmetric way by the terms — ^Vp a ln/„(r e ,p e ) and 

— |Vp- o ln/ y (r e , p e ), so that a numerical solution of these equations would be manifestly less computa- 
tionally intensive. 

Furthermore, because the equations ( pq ) are exact, they would be more reliable in ab initio calculations. 

Approximate solutions that arc valid through second order, may be achieved in a direct manner following 
the techniques outlined in p2| and may be written in the form of distorted and coupled Coulomb waves. 
e.g., 

U(fu, ?») = c nKt (Of.) = r(i - im-y-^xFi {imtX -»C*s) . (7i) 

where the conventional, two-body normalization procedure has been employed so that that incoming wave 
has unit magnitude. 

The complete three-body wavefunction may then be written in the form: 

*% C l C <fM = e^-f-.+a.-SOc {Ck a )C Vqa (C,JC_^ (C 9 J^ K7 (Ckj)C„ p {M (72) 

a product of five, kincmatically coupled Coulomb waves. To see that this solution is valid, first note that 
one can show in a straightforward manner that the distorted Coulomb waves have the following asymptotic 
form: 

C^G^-O^aj (73) 

Hence, due to the relationship (f70|), the asymptotic form in fio will be identical to that of the Redmond, 3C 
wavefunction. 



Note that the general form of the local distortion given in ( p2b| ) has been used to derive this form. 

1G 



The Jacobi Coordinates 




Figure 1: The Jacobi coordinates are most often used to study many-body kinematics, because any orthogonal 
pair. {f v , p v ), v = a, (3, 7 may be used. See equations (nf) for the definition of the reduced masses. 
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£ Dependence of the Local Distortion 
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Figure 2: The real part of the local distortion experienced by an electron-proton continuum pair with relative 
energies as indicated and a scattering angle of 9 — 4° . 
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H and Z Dependence of the Local Distortion 
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Figure 3: The real part of the local distortion experienced by both electron-electron and electron-proton 
continuum pairs, with values as in Fig. g. Note that the sign of the distortion is different for the attractive 
(e~,p + ) and repulsive (e~,e~) cases. 



Stationary Points of the Local Distortion 




Figure 4: The local distortion experienced by an electron-proton continuum pair with a relative energy of 
10 eV, and a scattering angle of 6 = 16°. Note that when the zeroes of real and imaginary part of the 
distortion coincide, the distortion contribution will be identically zero. 
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Topology of the Local Distortion 
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Figure 5: Shown here are (a) three-dimensional and (b) contour plots of the real part of the local distortion 
experienced by an electron-proton continuum pair with a relative energy of 10 cV. 
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Local Distortion Phase Effects 
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Figure 6: Shown here is the position dependent-phase, S Vp (^ p ). Note that it is very nearly constant over one 
log-cycle and that it is nearly independent of the relative energy of the electron-proton continuum pair. 



Local Distortion Effects for Various Continuum Pairs 
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Figure 7: Shown are the local distortion effects experienced by (a) electron-positron, (b) electron-electron, 
(c) electron-antiproton , (d) positron-proton, (e) proton-antiproton, and (f) proton-proton, continuum pairs 
with £ = 10 eV and 6 = 16°. While these effects become negligible asymptotically, the effects become more 
pronounced in the interior regions; the range has been shortened in (e) and (f) to show the dramatic change 
due to the increase in the reduced mass of the pair 
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WYJ Dalitz Plot £ t = 6.5 cV, V\ ov , = J% 




\ 

m 



(a) 



WYJ Dalitz Plot £ t = 7.5 eV, 7\ igh = J$ 
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Figure 8: These original data, obtained directly from the authors, show that as the total center-of-mass 
energy of the three-body system, £ tl is increased, a proportionately smaller number of particles undergoes 
kinematic rearrangement. Shown also are the predicted probability ratios for this rearrangement for the (a) 
£t = 6.5 eV and (b) £t = 7.5 eV, triple-coincidence events. 
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Explanation of Three-Body Distortion Effects 
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Coulomb Saddle 



Figure 9: (a) During the breakup of the H$ ion the HJ is emitted in the forward direction and the H%* is 
emitted in the backward direction, (b) The subsequent breakup of this ion occurs such that the H f and H~ 
ions are formed, and the H~ is preferentially associated with H f . (c) The continuum pair, (H t , H~) acquires 

a local momentum, P(f). pushing the H~ over the Coulomb Saddle, so that the H~ will be associated with 
the H t . (d) Because the two protons are in fact indistinguishable, one will observe an asymmetry between 
"low" and "high" energy scattering events. 
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Distortion Effects During (H 3 )* Breakup 
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Figure 10: During the breakup of the H$ ion, the magnitude of the local distortion effects experienced by 
the {Hf, H~) continuum pair while in the reaction zone are many orders of magnitude larger than the effects 
experienced by either the (e, e) or (e,e) continuum pairs. Here 6 = 176° and £ = 1.5 cV for comparison 
with the experiment. 
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Three-Body Gauge Variation with £ 
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Figure 11: The effects of the three-body gauge transformation are shown to depend critically upon the 
relative energy of the continuum pair. Observe that as the relative energy of an electron-proton continuum 
pair, with a scattering angle of 9 = 0°, is increased from (a) £ = 3 eV to (b) £ = 10 eV to (c) £ = 54.4 eV, 
the effects in the reaction zone nearly vanish. 
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Three-Body Gauge Variation with 9 
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Figure 12: The effects of the three-body gauge transformation are shown to depend critically upon the 
scattering angle of the continuum pair. Observe that as the scattering angle of an electron- proton continuum 
pair, with a relative energy of £ = 10 eV, is increased from (a) 9 = 0° to (b) 9 = 75° to (c) 9 — 150°, the 
effects in the reaction zone become more pronounced. 
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Three-Body Gauge Variation with /j, 
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Figure 13: The effects of the three-body gauge transformation are shown to depend critically upon the 
reduced mass of the continuum pair. Observe that as the reduced mass of the continuum pair, with a 
relative energy of £ = 10 cV and a scattering angle of 8 = 0°, is increased from (a) fi = \[i ev to (b) 
jjl = 10fi ep to (c) /i = 20/iep, the effects in the reaction are dramatically altered. 
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Asymptotic Phase Variation 
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Figure 14: Shown is the variation of the phase achieved by the three-body gauge transformation (see text 
equation ([43)) for (a) a proton-proton or proton-anti-proton, (b) an electron-proton or positron-proton and 
(c) electron-electron or electron-positron continuum pair. Observe that as the reduced mass of the system 
is decreased accordingly, the corresponding energy decrease is sufficient to detect a variation of one part in 
a million. 
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